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Abstract 

Classical hard spheres crystallize at equilibrium at high enough 
density. Crystals made up of stackings of 2-dimensional hexagonal 
close-packed layers (e.g. fee , hep , etc.) differ in entropy by only 
about W~ 3 kB per sphere (all configurations are degenerate in en- 
ergy). To readily resolve and study these small entropy differences, 
we have implemented two different multicanonical Monte Carlo algo- 
rithms that allow direct equilibration between crystals with different 
stacking sequences. Recent work had demonstrated that the fee stack- 
ing has higher entropy than the hep stacking. We have studied other 
stackings to demonstrate that the fee stacking does indeed have the 
highest entropy of all possible stackings. The entropic interactions 
we could detect involve three, four and (although with less statistical 
certainty) five consecutive layers of spheres. These interlayer entropic 
interactions fall off in strength with increasing distance, as expected; 
this fall-off appears to be much slower near the melting density than 
at the maximum (close-packing) density. At maximum density the en- 
tropy difference between fee and hep stackings is 0.001f5±0. 00004 ks 
per sphere, which is roughly 30% higher than the same quantity mea- 
sured near the melting transition. 
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1 Introduction 



The question of which crystalline stacking of hard spheres near close packing 
is thermodynamically stable is a long standing one. The interest is partly 
due to the extreme anharmonicity of hard core interactions and partly due 
to the fcc-hcp phase transition in solid helium HJ. This problem is diffi- 
cult both experimentally and theoretically. Experimentally, classical hard 
spheres are approximated by spherical colloidal particles with interactions 
whose ranges are very short compared to their radius. Deviations from ideal 
hard spheres are due to polydispersity of the spheres, and due to interactions. 
The van der Waals interaction can be reduced by matching the dielectric co- 
efficients of the particles and the solvent [|3], [| . Since, for ideal hard spheres, 
the free energy differences between the different stackings are very small, as 
we will see, one would expect the equilibration time to be very long. Most 
studies have seen a random stacking pattern. However, some experiments 
have reported that the observed random stacking patterns in slowly-grown 
or well-annealed colloidal crystals are biassed more towards fee rather than 
hep stacking 0. 

The free energy difference between different classical hard sphere crystals 
at fixed volume is only due to the entropy difference, since the energy is 
the same for all allowed configurations. The numerical work, before the 
present paper, had only looked at the hep and fee stackings. The first 
studies calculated pressure using molecular dynamics simulations and then 
obtained entropies by integrating the pressure vs. volume along reversible 
paths from states with known entropy || [7|, |S|, [5|. These studies were not 
able to detect the entropy difference between fee and hep crystals. Later, 



Frenkel and Ladd [ID instead integrated along a path connecting the hard 
sphere model to Einstein crystals of the same lattice structure, by adding 
to the model ideal springs tethering each ball to its lattice site. In this 
approach they integrate the derivative of free energy with respect to the 
spring constant. They came up with the bounds on the entropy difference 
per sphere: —0.001 < As* < 0.002 in units of k B , where As* = Sf cc — Sh cp - 
Recently, Bolhuis, Frenkel and the present authors used both a new im- 



plementation of the multicanonical Monte Carlo (MCMC) method \TA, 15 



as well as the Einstein crystal method to reduce the statistical errors down to 
the 10 _4 /cb per sphere level. This allowed us to accurately resolve the entropy 
difference of roughly 10 -3 fc# per sphere between fee and hep crystals, with 
the fee crystal having the higher entropy. [F2| This quantitatively corrected 
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the recent pressure-integration study of Woodcock JO] , confirming his result 
that the fee crystal has higher entropy. More recently, Bruce et al. have 
found a superior implementation of the multicanonical method for this prob- 
lem, reducing the statistical error in As* down to near the 10~ 5 /cb level |13| 



The various results for As* are summarized in Table |l[ It is clear from the 
table that MCMC is able to obtain substantially smaller statistical errors for 
this problem, compared to the more conventional integration methods. 

The reason why integration methods were used initially was that it did 
not appear possible for the hard-sphere system to equilibrate between the 
fee and hep crystal structures, due to the very large, or even infinite, free 
energy barrier separating them. This was certainly true for standard molec- 
ular dynamics or Monte Carlo methods. However, the MCMC method, the 
implementations of which we will summarize below, is designed precisely to 
eliminate such large free energy barriers and allow equilibration between very 
different states. This permits a direct measurement of the relative entropies 
of the two states simply by comparing the probabilities of their occurrences 
in a single simulation. 

Since only the hep and fee crystals had been examined in previous work, 
we have also looked at other stackings of hexagonally close packed planes 
of spheres to make certain that neglecting the other possible stackings was 
reasonable. The entropy differences between the stackings can be described 
as due to interactions between layers. We have been able to detect the 
entropic interaction between a given layer and its nearby layers that are two, 
three, and possibly even four layers away. These interactions are all of the 
sign that favors the fee stacking, so we confirm, as is no surprise, that the 
fee stacking has higher entropy than all other stackings, not just higher than 
hep . At the maximum packing density, we find that the interaction with the 
third-neighbor layer is roughly an order of magnitude smaller than that with 
the second-neighbor layer, as seems quite reasonable. For lower density, near 
the melting transition, the fall-off of these entropic interactions with distance 
appears to be much slower, presumably reflecting the larger fluctuations of 
the individual sphere positions. 

For any given stacking, the entropy varies as a function of any homoge- 
neous lattice deformation at constant volume fraction. For the fee stacking, 
the undeformed lattice has cubic symmetry, so must by symmetry be a sta- 
tionary point of the entropy vs. deformation, and it is the maximum. For 
the other stackings, there is no such symmetry, and the maximum entropy 
may be obtained for a deformation where the expansion of the lattice away 
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from close-packing is not isotropic. We have looked for this possible effect in 
the hep stacking by measuring the entropy vs. the uniaxial lattice deforma- 
tion (the c/a ratio). If there is an anisotropy, we were unable to detect it. 
If there is an entropy difference between the highest entropy state and the 
isotropically-expanded state for the hep stacking, this difference is less than 
10~ 5 ks per sphere, so is well below the statistical errors in our comparisons 
to the other stackings. 

Another issue that arises in simulating hard-sphere crystals is whether 
collisions between spheres that are not nearest-neighbors can be neglected. 
It is certainly more convenient to make the approximation of including the 
hard-sphere interaction only between nearest-neighbors. Of course this ap- 
proximation is terrible in the liquid phase, but we have found that it is 
actually quite good in the solid phase even at the melting density. There we 
detected no difference in As* between the model where only nearest-neighbor 
spheres interact and the model where second-neighbors also interact, indi- 
cating that the difference is also smaller than our statistical errors when 
comparing the entropies of different stackings. It is usually not clear in the 
literature whether or not further neighbor interactions have been taken into 
account. This result of ours shows that it doesn't matter within the solid 
phase at the presently available resolution of the entropy. 

2 Model 

The model we study is hard spheres: classical monodisperse spheres that 
are forbidden to overlap. All permitted configurations (with no overlaps 
between spheres) have the same energy, which we may set as zero energy. 
At high enough density this system crystallizes at equilibrium, and it is this 
crystalline equilibrium phase that we study here. Some of our results are 
for the maximum possible, or close-packed density, which means the system 
is being treated perturbatively to lowest order in the difference between the 
density and the close-packed density. In this limit the system is equivalent 
to a simpler system of aligned, hard dodecahedra 0. 

We consider close-packed crystal structures that consist of planes of 
hexagonally close-packed (in 2-dimensions) spheres stacked up in the ver- 
tical direction. As is standard in discussing close-packed crystals, the stack- 
ing sequence can be denoted by a sequence of the letters A, B and C, with 
nearest-neighbor layers in the sequence always having different letters. Any 



4 



global permutation on the letters in the sequence simply represents a rotation, 
reflection or translation of the structure, so will not change the entropy. AB- 
CABCA... is a sequence that represents the fee stacking, while ABABABA... 
represents the hep stacking. To fully remove the degeneracy associated with 
the permutations, we may assign a Ising-like spin Oi to each layer % based on 
the local stacking sequence of that layer and its nearest-neighbor layers im- 
mediately above and below it. If that local stacking matches the fee pattern 
(e.g., ABC), <7j = +1, while Oi = —1 if it instead matches the hep pattern 
(e.g., ABA). For the three layers being compared, all local stacking pat- 
terns are equivalent to one of these two under permutations of the labels 
A, B and C. This maps each stacking sequence on to a spin pattern of a 
one-dimensional Ising model. 

The entropy of a given stacking is a function of the stacking sequence, 
which is described by the spins Uj. It is reasonable to expect that the shortest- 
range entropic interactions are the largest, so the total entropy may be ex- 
panded as 

n n n n 

S = Ans + Ah (Ti + AJ + AJ' ^ + Ah' ^ <JiOi+\<Ji+i H , 

i=l i=l i=l i=l 

(1) 

for a stack of n layers containing A spheres per layer. We have periodic 
boundary conditions so <7 n +i = a±, etc. The bulk of the entropy is inde- 
pendent of the stacking sequence and given by s per sphere; s is strongly 
density-dependent. The shortest-range entropic interaction is h, which in- 
volves the sequence over three consecutive layers of spheres; this is the short- 
est sequence that can have distinct stackings. This term is the magnetic field 
in the corresponding one-dimensional Ising model. The next term is the in- 
teraction J between adjacent spins, and arises from the entropic interactions 
among four consecutive layers of spheres that are not already captured by 
the first term h. We find, as is reasonable, that J < h. The next longer-range 
interactions {h! and J 1 ) that involve five consecutive layers are also displayed 
above; only in one case (h' for density near melting) could we detect these 
interactions in our simulations at a level that may be statistically significant. 

The intuition behind this model is that the entropy of a sphere is mostly 
determined by how it is caged by its nearest neighbors and to a progres- 
sively lesser extent by the further neighbors. The interaction parameters do 
depend on the density. We find that, in units of 10 _5 /cs per sphere, the en- 
tropic interactions that we could detect change from (h, J) = (55, 6) at the 
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highest density (close-packing) to (h, J, h!) = (37, 18, 9) at the lowest density 
that the equilibrium crystal can have before it melts (at roughly 74% of the 
close-packed density) . All detected interactions are of the sign such that the 
fee stacking has the largest entropy. 



3 Multicanonical Monte Carlo Method 

To make direct comparisons of the entropies of hard-sphere crystals with 
different stacking sequences, we want an algorithm that will produce direct 
equilibration between the two sequences. Then the entropy difference is 
simply the logarithm of the ratio of the equilibrium probabilities of the system 
exhibiting the two sequences in question. Of course, near close-packing the 
hard spheres in a physically realistic molecular dynamics or Monte Carlo 
simulation are strongly trapped by their neighbors, so the stacking pattern 
will not change in any reasonable time scale. The multicanonical method 
b|, [L5| was invented to allow systems to transform at equilibrium between 



states that are separated by a high free energy barrier. This method has 
been generalized and applied to this hard-sphere crystal problem in two ways, 
which we describe next. 

In both implementations the position of sphere i is given as n = Ri + Ui, 
where Ri is the ideal reference lattice position in the absence of fluctuations 
and Ui is the displacement of sphere i away from that reference position. The 
algorithms both have Monte Carlo moves that change the reference lattice 
without changing the displacements Ui, as well as more standard moves that 
move individual spheres without changing the reference lattice. We describe 
the shear implementation first; this method we developed and used to obtain 
our first results flT^f . However, Bruce, et al. jn| subsequently developed 
the simpler overlap implementation, which we find is computationally more 
efficient and easier to program, so we used it for all of our more recent 
simulations. 



3.1 Shear Implementation 

In the shear implementation we used an equally-spaced sequence of ideal 
reference lattices, Ri(A), labelled by an index, A = 0, 1, 2, h, that linearly 
interpolate between the two stackings of interest, which are the beginning 



6 



(A = 0) and the end (A = h) of that sequence. Thus we have 

hRi(X) — (h — A)Ri(O) + XRi(h). (2) 

This produces a A-dependent relative shear between each pair of adjacent 
layers whose local stacking pattern changes. The reference sites in the inter- 
mediate lattices (0 < A < h) are not all equally spaced, and generally some 
are too close together for the hard spheres to fit without touching. In the 
original model, two spheres are assumed to touch if their separation, |n — rj|, 
is less than d, the diameter of a sphere, and this must remain true for the 
two stackings of interest, which are represented by A = and h. For the 
intermediate lattices, on the other hand, we allow the distance of contact, 
dij(X) to be either larger or smaller for each pair of nearby spheres in adjacent 
layers whose relative reference position, Ri — Rj, changes with changing A. 
We attempt, using feedback, to choose these dy(A) so that the entropy is a 
monotonic function of A and the average displacements, < Ui > vanish for 
all i and A. Note that the pairwise interactions are different for each interpo- 



lation point, A. This is different from the original MCMC method |14], [L5|, 
where the Gibbs distribution is multiplied by a A-dependent (but otherwise 
configuration-independent) reweighting factor, in order to make the free en- 
ergy monotonic between the two states of interest, thus eliminating the free 
energy barrier. 

To start simulating one has to choose how many interpolation points to 
use, h, and values for the dij(X). There are two types of moves. One is a 
single sphere move, changing one of the displacements Ui. The other is a A- 
move that increases or decreases A by one without changing any of the sphere 
displacements. Any attempted move of either type is accepted if it does not 
result in any contact between spheres. The entropies and average displace- 
ments, as well as the acceptance rates of the moves are measured. Based on 
these results h is adjusted to attempt to maximize the rate of equilibration 
between the two stackings of interest, and the rfy(A) are adjusted to attempt 
to make the entropy monotonic and eliminate the average displacements. If 
this feedback is successful, we then measure the relative entropy of the two 
stackings. 

We succeeded in getting this procedure to work for comparing the fee and 
hep stackings for lattices of size up to 8 3 . However, the difficulty of getting 
the feedback to converge appeared to be increasing strongly with lattice 
size. Typically, "bottlenecks" would form between the two stackings of in- 
terest where the A-move acceptance rate was very small or zero, preventing 
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equilibration, and the attempts at eliminating these bottlenecks through the 
feedback were time-consuming and not always succesful. However, we were 
able to obtain the entropy difference between the fee and hep stackings to 
within statistical errors of roughly 10 _4 fc# per sphere [ 12f , as is summarized 
in Table 1. Then we learned of the much more staightforward overlap imple- 



mentation of MCMC for this problem reported by Bruce, et al. 
we discuss next. 



HJ, which 



3.2 Overlap Implementation 



The overlap implementation of MCMC [13[ uses only two reference lat- 
tices, which are the two different stackings whose entropy we are comparing. 
Let us call these two reference lattices a and (3. Again there are standard 
single-sphere moves and changes of the reference lattice. For any reasonable 
sized lattice, the move that changes reference lattices will be rejected due to 
sphere overlaps for all but an infinitesimal fraction of the sphere configura- 
tions. What is needed is to bias the simulation towards those rare sphere 
configurations that allow the stacking to be changed. 

To do this, Bruce et al. [O introduce the overlap order parameter 



M{{u}) = M({u}, a) - M({u},(3) (3) 

where M({w},7) is the number of pairs of spheres that overlap in the con- 
figuration {u} for stacking 7. For any allowed configuration with stack- 
ing 7, M({m},7) = 0, but for configurations of the other stacking, usually 
M({u}, 7) > 0. To have the change of stacking be an allowed move, we need 
({«}-) = 0, so no overlaps are produced by the move that changes the 
reference lattice. 

The overlap multicanonical simulation samples the biassed, but unnor- 
malized, distribution 

P({u} n \{ V }) = P({u} n )e^ M ^\ (4) 

where P({w},7) for (unbiassed) hard spheres is simply a constant for all al- 
lowed sphere configurations in stacking 7, and is zero otherwise. The weights 
{77} are chosen to eliminate the free energy barrier separating the two stack- 
ings, thus allowing equilibration between them. 

Let P(Ai) be the equilibrium, normalized, unbiassed (t)(M) = for all 
M) probability distribution of the overlap, assuming the system does fully 
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equilibrate between the two stackings. Then the probability of being in 
stacking a is 



P 



P(0) 



+ E P{M) 



(5) 



a 



2 



M<0 



and the entropy difference we are interested in is 



kg In 



( 



P 



a 



) 



(6) 



l-P 



a 



P(M) has a local minimum at A4 = and local maxima at positive and 
negative A4 that represent the most probable overlaps for the two differ- 
ent stackings. The weights f](A4) are chosen to be nonzero only between 
these two local maxima of P(M). The unnormalized, biassed probability 
distribution for A4 is 



We choose the weights {rj} so that P(A4 \ {rj}) is linear in A4 between 
the maxima of P(A4). A simulation with a given set of weights produces 
estimates of P(A4) and also a new estimate of what the appropriate weights 
are to achieve this linearity. These new weights are then used for the next, 
longer simulation if the statistical errors have not yet been reduced down to 
the desired level. This procedure straightforwardly and effectively eliminates 
the free energy barrier between the two stackings and allows an accurate 
measurement of the entropy difference. 

The overlap implementation of MCMC does not suffer from the ten- 
dency to form "bottlenecks" that slowed down the equilibration between 
the two stackings in our shear implementation of MCMC. We used the over- 
lap method to obtain most of the data reported in this paper. Where we 
compared the two implementations the measured entropy differences were, 
of course, the same. 

3.3 Boundary Conditions 

Suppose one stacks N 3 planes of 2-dimensionally (hexagonally) close-packed 
spheres to form an arbitrary stacking. Each plane has iVi x N2 spheres Q, 

^^We always choose Ni — N2 to preserve the hexagonal in-plane symmetry 



P(M I M) 



P{M) e^ M \ 



(7) 
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Ni in the i-direction and A 2 in the j-direction. The z-direction is chosen to 
coincide the the x-direction and the j-direction is chosen at 60° anti-clockwise 
from the x-direction looking from above, i and j are the two basis directions 
of the 2-dimensional close packing. Each site has coordinate (i, j) in-plane, 
z = 0, 1, . . . , ^ - 1 and j = 0, 1, . . . , N 2 - 1. 

Any stacking can be formed by fixing the position of the k) reference 
site (in layer k) relative to the same site (i,j,k ± 1) in the nearest layers. 
Define S fc = R(i,j,k + 1) — R(i,j,k), where R are the reference sites. We 
take our unit length to be the lattice spacing so = 1. We use 

if from layer k to layer k + 1 is a 

forward permutation of ABC, , , 

if from layer k to layer k + 1 is a 
backward permutation of ABC. 

Taking = S + for all k gives an fee stacking. For the hep stacking £& = S + 
for k even and = XL for k odd (or vice versa). 

All simulations are done with periodic BC which is implemented in the 
usual way: 

sphere at site (i, j, k) = sphere at site (i + Ni,j, k) 

= sphere at site + N 2 , k) 
= sphere at site k + N 3 ), 

so S fc = Sfc + 7v 3 . Our implementation of periodic BC allows any A3 which is 
a multiple of both of the periods of the patterns of the two stackings in a 
given simulation. For example, the period for fee is one layer, while that for 
hep is two layers, so A3 can be any even number when we compare fee and 
hep entropies. 

4 Results 

Our results for entropy differences between different stackings at both close- 
packing (p/p C p = 1) and near melting (p/p cp = 0.739) for different system 
sizes are summarized in Table 0. A statistically significant finite-size effect 
was detected only for the smallest size 4 3 . Based on this, we assume that the 
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finite-size effect is negligible for sizes 8 3 and larger, and we use those data 
for calculating the values of the entropic interactions using our Eq.[L| 

At close packing we fit, using the results of N = 8 3 , 9 3 , and 10 3 and the 
various different stackings, first letting all four parameters {h, J, J',h'} vary, 
and then setting h! = J' = and only varying {h, J}. In the former case 
we obtain h = 54.6 ± 2.8, J = 6.1 ± 1.6, J' = -0.3 ± 1.1 and h' = 3.4 ± 2.7 
in units of lO" 5 ^ with x 2 /(df = 3) = 1.5 and for the ti = J' = fits: 
h = 57.2 ± 2.1 and J = 6.0 ± 2.2 with x 2 /(df = 5) = 1.9. We can see that 
the first fit with J' and h! allowed to vary gives values of them consistent 
with zero. Comparing the two fits we also see that the inclusion of these 
two longer-range interactions in the fit does not significantly perturb the 
values of the shorter-range interactions h and J. We therefore conclude 
that our data can be explained using the model with only h and J non-zero. 
Indeed, J' and h! not being important is consistent with the small system size 
N = 8 3 being close to the thermodynamic limit. The signs of the non-zero 
interactions all favor fee stacking, therefore fee has higher entropy than all 
other stackings, consistent with experiment 0]. Notice that h » J » J' , h' 
shows that the entropic interactions decrease rapidly as their range increases. 

Similar fits of our data at a density (p/p cp = 0.739) near melting yield 
h = 36.9 ± 3.1, J = 18.2 ± 3.0, J' = 2.5 ± 2.2 and h! = 8.8 ± 2.8 with 
X 2 / {df = 2) = 0.2. Thus it appears that the entropic interactions decrease in 
relative magnitude with distance much more slowly at this lower density than 
they do near close-packing, which is perhaps expected, given the larger free 
volume which allows the spheres to make larger excursions away from their 
ideal lattice positions. Again, the interactions are all of the sign that favor 
the fee stacking. Our detection of h! is only at the three standard deviation 
level, so has a small chance of being just a statistical fluctuation in the data. 
However, if we fit assuming h' — J' — the quality of the fit declines strongly, 
giving h = 43.2 ± 3.7 and J = 17.3 ± 3.8 with x 2 /(df = 4) = 5.4. 

For a general stacking pattern, the expansion as the density is reduced 
from close-packing need not be isotropic. For the fee stacking it must, by 
the cubic symmetry, but for the other stackings, the expansion along the 
direction normal to the layers can be different from that along the directions 
parallel to the layers. We have tested for this by allowing the ratio of these 
two expansions to vary in a simple simulation of the hep stacking at close- 
packing, measuring the entropy vs. the ratio, and fitting to find the ratio 
that maximizes the entropy. We find that this optimal ratio is within ±0.002 
of isotropic, and the entropy difference between isotropic expansion and the 
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optimal expansion ratio is no more than 10~ 5 ks per sphere, so is smaller 
than the statistical uncertainties in our simulations. Because of this we have 
always assumed isotropic expansion in the the entropy comparisons we have 
made. 

The issue of further neighbor interactions arises for p/p cp < 1. At close- 
packing it suffices to test only for collisions between nearest-neighbor spheres 
because further neighbors cannot touch. Not testing for further-neighbor 
collisions speeds up the computer program. As the density is reduced can we 
keep this approximation? For a model that includes only nearest-neighbor 
interactions between spheres, the crystal is actually only metastable: once 
a sphere "escapes" from the cage of its nearest-neighbors it wanders freely. 
We find that for densities at or above the melting density, the rate of these 
"escapes" is very low, allowing a good measurement of the entropy of the (now 
metastable) crystal. We have also measured the entropy differences between 
the model with only nearest-neighbor interactions and the otherwise identical 
model with nearest and next-nearest-neighbor interactions, S n — S nn , near 
melting (see Table 0). Of course, adding the extra interactions does reduce 
the entropy a little (roughly 8 x 10~ 5 ks per sphere), but this reduction is the 
same, within errors, for both fee and hep stackings. Thus we conclude that 
any systematic error in our entropy comparisons due to using only nearest- 
neighbor interactions are smaller than the statistical errors. Therefore, we 
have used the faster nearest-neighbor only model in most of our simulations 
near the melting density. 
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MCMC, overlap implementation 
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MCMC, overlap implementation 


present 


1.00 
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(20) 


MCMC, shear implementation 


present 


0.731 


512 


85 


(10) 


MCMC, shear implementation 


present 



Table 1: Recent results of fee- hep simulations for various densities (scaled 
by the close-packed density p cp ). N is the number of spheres in the samples 
simulated. The entropy difference per sphere is As*, with the statistical 
errors in parenthesis, (fee has higher entropy.) Please note that the errors 
are particularly small for the overlap implementation of MCMC developed 
by Bruce, et al. 
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61.2 


(2.2) 


1 


8 3 




— s_ 






h — 2J + h' 


44 


(4) 


1 


8 3 s +++ 


+++ — 




h-++- 




\J+\J' + h' 


8.2 


(2.2) 


0.739 


8 3 


s+ 


— s_ 






2h + 2h! 


90.2 


(4.3) 


0.739 


8 3 




- - s+- 






J-2J' 


13.7 


(2.9) 


0.739 


8 3 


s+- 


— s_ 






h — 2J + h' 


10.5 


(5.0) 


0.739 


8 3 s +++ 


+++ — 




h-++- 




\j+\j' + ti 


19 


(3) 


0.739 


8 3 s+ 


++++- 


f- - s+- 


h+- 




\J> + \h> 


6.2 


(3.2) 


0.739 


8 3 s +++ 


+-++- 


- 8++— 


--++- 




\h + \J>-\h> 


18.0 


(1.9) 



Table 2: This table summarizes the entropy differences per sphere between 
various pairs of stackings at the densities we studied. The subscript on s is a 
sequence of cr^'s that is repeated to get the stacking sequence, so + denotes 
fee , — denotes hep , and the others are less simple stackings (see text). 
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P/Pcp 


N As 


10 5 


x As/k B 


0.739 


oiS fee _ fee 
u n J nn 


8.3 


(1.9) 


0.739 


o3 Q hcp Q hcp 
° b n b nn 


7.8 


(1.9) 



Table 3: The entropy differences per sphere between hard sphere crystals with 
only nearest neighbor interactions (n) and the otherwise identical system with 
both nearest and next nearest neighbor interactions (nn) . The systems with 
added interactions have lower entropy, but the change is independent of the 
stacking pattern at the present statistical accuracy. 
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